Molecular and Clinical Epidemiology of SARS-CoV-2 Infection Among Vaccinated and Unvaccinated Individuals in a Large Healthcare Organization from New Jersey (preprint)

New Jersey was among the first states impacted by the COVID-19 pandemic, with one of the highest overall death rates in the nation. Nevertheless, relatively few reports have been published focusing specifically on New Jersey. Here we report on molecular, clinical, and epidemiologic observations from the largest healthcare network in the state, in a cohort of vaccinated and unvaccinated individuals with laboratory-confirmed SARS-CoV-2 infection. We conducted molecular surveillance of SARS-CoV-2-positive nasopharyngeal swabs collected in nine hospitals from December 2020 through June 2022, using both whole genome sequencing (WGS) and a real-time RT-PCR screening assay targeting spike protein mutations found in variants of concern (VOC) within our region. De-identified clinical data were obtained retrospectively, including demographics, COVID-19 vaccination status, ICU admission, ventilator support, mortality, and medical history. Statistical analyses were performed to identify associations between SARS-CoV-2 variants, vaccination status, clinical outcomes, and medical risk factors. A total of 5,007 SARS-CoV-2-positive nasopharyngeal swabs were successfully screened and/or sequenced. Variant screening identified three predominant VOC, including Alpha (n =714), Delta (n =1,877), and Omicron (n =1,802). Omicron isolates were further sub-typed as BA.1 (n =899), BA.2 (n =853), and BA.4/BA.5 (n =50); the remaining 614 isolates were classified as “Other”. Approximately 31.5% (1,577/5,007) of the samples were associated with vaccine breakthrough infections, which increased in frequency following the emergence of Delta and Omicron. Severe clinical outcomes included ICU admission (336/5007 = 6.7%), ventilator support (236/5007 = 4.7%), and mortality (430/5007 = 8.6%), with increasing age being the most significant contributor to each (p <0.001). Unvaccinated individuals accounted for 79.7% (268/336) of ICU admissions, 78.3% (185/236) of ventilator cases, and 74.4% (320/430) of deaths. Highly significant (p <0.001) increases in mortality were observed in individuals with cardiovascular disease, hypertension, cancer, diabetes, and hyperlipidemia, but not with obesity, thyroid disease, or respiratory disease. Significant differences (p <0.001) in clinical outcomes were also noted between SARS-CoV-2 variants, including Delta, Omicron BA.1, and Omicron BA.2. Vaccination was associated with significantly improved clinical outcomes in our study, despite an increase in breakthrough infections associated with waning immunity, greater antigenic variability, or both. Underlying comorbidities contributed significantly to mortality in both vaccinated and unvaccinated individuals, with increasing risk based on the total number of comorbidities. Real-time RT-PCR-based screening facilitated timely identification of predominant variants using a minimal number of spike protein mutations, with faster turnaround time and reduced cost compared to WGS. Continued evolution of SARS-CoV-2 variants will likely require ongoing surveillance for new VOCs, with real-time assessment of clinical impact.

Rare predicted loss-of-function variants of type I IFN immunity genes are associated with critical COVID-19 (preprint)

Background We previously reported inborn errors of TLR3- and TLR7-dependent type I interferon (IFN) immunity in 1-5% of unvaccinated patients with life-threatening COVID-19, and autoantibodies against type I IFN in another 15-20% of cases. Methods We report here a genome-wide rare variant burden association analysis in 3,269 unvaccinated patients with life-threatening COVID-19 (1,301 previously reported and 1,968 new patients), and 1,373 unvaccinated SARS-CoV-2-infected individuals without pneumonia. A quarter of the patients tested had antibodies against type I IFN (234 of 928) and were excluded from the analysis. Results No gene reached genome-wide significance. Under a recessive model, the most significant gene with at-risk variants was TLR7, with an OR of 27.68 (95%CI:1.5-528.7, P=1.1x10-4), in analyses restricted to biochemically loss-of-function (bLOF) variants. We replicated the enrichment in rare predicted LOF (pLOF) variants at 13 influenza susceptibility loci involved in TLR3-dependent type I IFN immunity (OR=3.70 [95%CI:1.3-8.2], P=2.1x10-4). Adding the recently reported TYK2 COVID-19 locus strengthened this enrichment, particularly under a recessive model (OR=19.65 [95%CI:2.1-2635.4]; P=3.4x10-3). When these 14 loci and TLR7 were considered, all individuals hemizygous (n=20) or homozygous (n=5) for pLOF or bLOF variants were patients (OR=39.19 [95%CI:5.2-5037.0], P=4.7x10-7), who also showed an enrichment in heterozygous variants (OR=2.36 [95%CI:1.0-5.9], P=0.02). Finally, the patients 13 with pLOF or bLOF variants at these 15 loci were significantly younger (mean age [SD]=43.3 [20.3] years) than the other patients (56.0 [17.3] years; P=1.68x10-5). Conclusions Rare variants of TLR3- and TLR7-dependent type I IFN immunity genes can underlie lifethreatening COVID-19, particularly with recessive inheritance, in patients under 60 years old.

Post-vaccination SARS-COV-2 among healthcare workers in New Jersey: a genomic epidemiological study (preprint)

ABSTRACT Emergence of SARS-CoV-2 with high transmission and immune evasion potential, the so-called Variants of Concern (VOC), is a major concern. We describe the early genomic epidemiology of SARS-CoV-2 recovered from vaccinated healthcare professionals (HCP). Our post-vaccination COVID-19 symptoms-based surveillance program among HCPs in a 17-hospital network, identified all vaccinated HCP who tested positive for COVID-19 after routine screening or after self-reporting. From 01/01/2021 to 04/30/2021, 23,687 HCP received either mRNA-1273 or BNT162b2 mRNA vaccine. All available post-vaccination SARS-CoV-2 samples and a random collection from non-vaccinated patients during the similar timeframe were subjected to VOC screening and whole genome sequencing (WGS). 62% (23,697/37,500) of HCPs received at least one vaccine dose, with 95% (22,458) fully vaccinated. We detected 138 (0.58%, 138/23,697) COVID-19 cases, 105 among partially vaccinated and 33 (0.15%, 33/22,458) among fully vaccinated. Five partially vaccinated required hospitalization, four with supplemental oxygen. VOC screening from 16 fully vaccinated HCPs identified 6 (38%) harboring N501Y and 1 (6%) with E484K polymorphisms; concurrent non-vaccinated samples was 37% (523/1404) and 20% (284/1394), respectively. There was an upward trend from January to April for E484K/Q (3% to 26%) and N501Y (1% to 49%). WGS analysis from vaccinated and non-vaccinated individuals indicated highly congruent phylogenies. We did not detect an increased frequency of any RBD/NTD polymorphism between groups (P>0.05). Our results support robust protection by vaccination, particularly among recipients of both doses. Despite VOCs accounting for over 40% of SARS-CoV-2 from fully vaccinated individuals, the genomic diversity appears to proportionally represent those among non-vaccinated populations. IMPORTANCE A number of highly effective vaccines have been developed and deployed to combat the COVID-19 pandemic. The emergence and epidemiological dominance of SARS-CoV-2 mutants, with high transmission potential and immune evasion properties, the so-called Variants of Concern (VOC), continues to be a major concern. Whether these VOCs alter the efficacy of the administered vaccines is of great concern, and a critical question to study. We describe the initial genomic epidemiology of SARS-CoV-2 recovered from vaccinated healthcare professionals and probe specifically for VOC enrichment. Our findings support the high-level of protection provided by full vaccination despite a steep increase in the prevalence of polymorphisms associated with increased transmission potential (N501Y) and immune evasion (E484K) in the non-vaccinated population. Thus, we do not find evidence of VOC enrichment among vaccinated groups. Overall, the genomic diversity of SARS-CoV-2 recovered post-vaccination appears to proportionally represent the observed viral diversity within the community.

Emergence of multiple SARS-CoV-2 antibody escape variants in an immunocompromised host undergoing convalescent plasma treatment (preprint)

SARS-CoV-2 Variants of Concerns (VOC), e.g., B.1.351 (20H/501Y.V2) and P1 (20J/501Y.V3), harboring N-terminal domain (NTD) or the receptor-binding domain (RBD) (e.g., E484K) mutations, exhibit reduced in vitro susceptibility to convalescent serum, commercial antibody cocktails, and vaccine neutralization, and have been associated with reinfection. The accumulation of these mutations could be the consequence of intra-host viral evolution due to prolonged infection in immunocompromised hosts. In this study, we document the microevolution of SARS-CoV-2 recovered from sequential tracheal aspirates from an immunosuppressed patient on tacrolimus, steroids and convalescent plasma therapy, and identify the emergence of multiple NTD and RBD mutations associated with reduced antibody neutralization as early as three weeks after infection. SARS-CoV-2 genomes from the first swab (Day 0) and three tracheal aspirates (Day 7, 21 and 27) were compared at the sequence level. We identified five different S protein mutations at the NTD or RBD regions from the second tracheal aspirate sample (21 Day). The S:Q493R substitution and S:243-244LA deletion had ~70% frequency, while ORF1a:A138T, S:141-144LGVY deletion, S:E484K and S:Q493K substitutions demonstrated ~30%, ~30%, ~20% and ~10% mutation frequency, respectively. However, the third tracheal aspirate sample collected one week later (Day 27) was predominated by the haplotype of ORF1a:A138T, S:141-144LGVY deletion and S:E484K (> 95% mutation frequency). Notably, S protein deletions (141-144LGVY and 243-244LA deletions in NTD region) and substitutions (Q493K/R and E484K in the RBD region) previously showed reduced susceptibly to monoclonal antibody or convalescent plasma. The observation supports the hypothesis that VOCs can independently arise and that immunocompromised patients on convalescent plasma therapy are potential breeding grounds for immune-escape mutants.